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Effects of vacuum polarization by hadronic and heavy-fermion insertions were the last unknown two-loop 
QED corrections to high-energy Bhabha scattering and have been first announced in [lj]. Here we describe the 
corrections in detail and explore their numerical influence. The hadronic contributions to the virtual 0(a 2 ) 
QED corrections to the Bhabha-scattering cross-section are evaluated using dispersion relations and computing 
the convolution of hadronic data with perturbatively calculated kernel functions. The technique of dispersion 
integrals is also employed to derive the virtual 0(a ) corrections generated by muon-, tau- and top-quark loops 
in the small electron-mass limit for arbitrary values of the internal-fermion masses. At a meson factory with 1 
GeV center-of-mass energy the complete effect of hadronic and heavy-fermion corrections amounts to less than 
0.5 per mille and reaches, at 10 GeV, up to about 2 per mille. At the Z resonance it amounts to 2.3 per mille 
at 3 degrees; overall, hadronic corrections are less than 4 per mille. For ILC energies (500 GeV or above), the 
combined effect of hadrons and heavy fermions becomes 6 per mille at 3 degrees; hadrons contribute less than 
20 per mille in the whole angular region. 

PACS numbers: ll.15.Bt, 12.20.Ds 
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I. INTRODUCTION 



Elastic e + e scattering, or Bhabha scattering, 



(Pi) 



(Pa) 



(P3) + e + (p 4 ) , 



(1) 



was one of the first scattering processes that were observed and predicted in quantum mechanics J2]. It has 
a unique and clean experimental signature. The accuracy of theoretical predictions profits from its purely 
leptonic external particle content and from the extremely small electron mass. The first complete one-loop 
prediction in the Standard Model was yD, the first Oia) predictions in the Standard Model with account of 
hard bremsstrahlung were determined in HHIaLZl&Eali me effects from hadronic vacuum polarization were 
first studied in 11011 . and the leading NNLO corrections from the top quark in lUlll . The complete electroweak 
two-loop corrections are available in form of few form factors 11121 [1311 . but they are not implemented for 
Bhabha scattering so far. During the years, a rich literature on the subject arose, both concerning QED Monte 
Carlo results and virtual electroweak corrections; see II 1 41, 1 1 5l 1 1 6l 1 1 7U 1 8l 1 1 9L |2.(A 12. lLl22l J23L J24l J25L )26f 
2§ [3^ 1J IIlS JllJIlJI^]^^ 



£8i i5j, i6oi |6H i62j, mm mm m m mMuMuIl MMm 

also the references therein. 

Quite recently, an experimental precision at the per mille level or beyond seems feasible both at meson 
factories and in the ILC (and GigaZ) project JH, ill IM HHS IH . As a reaction to that, a program of 
systematic evaluation of t he complete next-to-next-to leading order (NNLO) contributions was emerging ll 



IMRI<^IMI<rMI<^[M[ioil[^ 



In this article, we extensively describe the evaluation of the last building block of QED two-loop corrections, 
namely the corrections from heavy fermi ons and hadronic vacuum polarization. Note that the latter result 
has been confirmed very recently in lll40ll (upon using the same parametrisation of the vacuum polarization, 
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the agreement between the two studies is perfect, 5 digits for the 0{a A ) NNLO terms). Both for reasons of 
completeness and in order to ensure easy comparisons, we will also include in the discussion the Nf = 1 
corrections which consist of purely photonic corrections and electron loop insertions, the soft bremsstrahlung 
and soft electron pair emission corrections. All the two-loop contributions are calculated in our numerical 
Fortran package bhbhnnlohf . F and will be made available at the webpage 1 1 1*711 - 

The organization of the paper is as follows. In Section[n] we introduce notations and the Born cross-section. 
Section [Til] collects the known facts on pure vacuum-polarization corrections as they will be used, and Sec- 
tion|IV]the pure self-energy corrections to the cross-section. Section[V]contains the irreducible vertex correc- 
tions and Section [VT] the various infrared divergent corrections, including reducible corrections, soft-photon 
emission and the most complicated ones from the irreducible two-loop box diagrams. The three kernel func- 
tions for the latter have been evaluated for the first time. Section [VTIl contains a discussion of numerical effects 
at a variety of energies, typical of meson factories, LEP, ILC. In the Summary we will also point to potential 
further research. AppendiceslAltolFlare devoted to technical details of fermionic vacuum polarization, one-loop 
master integrals, soft real bremsstrahlung, real pair emission, the evaluation of the hadronic cross-section ratio 
i?had, and on our evaluation of complex polylogarithms. Some Mathematica files of potential public interest 
and the Fortran package are available at the webpage Ill7ll . 



II. THE BORN CROSS-SECTION 

The QED tree-level differential Bhabha-scattering cross section with respect to the solid angle il, in the 
kinematic region m 2 <C s, \t\, \u\, is: 



da a 2 fvi(s,t) | n V2(s,t) t Vi(t,s) 

V st + t 2 

2 / +\ 2 

as t x 



tU +1+ ;'- (2) 



Here, a is the fine-structure constant 111811 . 

a = 1/137.035999679(94), (3) 

and 

vi(x,y) = x 2 + 2y 2 + 2xy, (4) 

V2(x,y) = {x + y) 2 . (5) 

The cross-section depends on the Mandelstam invariants s, t and u, which are related to E, the incoming- 
particle energy in the center-of-mass frame, and 9, the scattering angle: 

s = (Pi+p 2 ) 2 = 4£ 2 , 

t = ( Pl -p 3 ) 2 =-<iE 2 S in 2 ^y 

u = ( Pl -p 4 f = -4E 2 cos 2 (J\ , (6) 

where 

s + t + u = 0. (7) 

For the numerical estimates at higher energies, it is reasonable to normalize the higher order corrections to 
the complete electroweak effective Born cross-section: 



da K 



a 2 



dn 

with 



— (T s + T st + T t ) , (8) 



T s = (l + cos 2 0) l + 2Re X (s)(v 2 ) + \x(s)\ 2 {l + v 2 ) 2 + 2cos0 [2Re X (s) + |x(s)| 2 M] , (9) 



(1 + COS60 2 
T * = ~ 2 71 ^r{l + [xW+Rex(.s)](l + w 2 )+x(i)Rex(s)[(l + « 2 ) 2 +4« 2 ]}, (10) 



(1 — cos( 
(l + cosf^ 2 
(1 — cos( 



r * = 2 ;; I Zl i 1 + 2 *w (i + . 2 ) + xw 2 [d + v-f + 4, 2 ] } 



(l-cos6») 2 
We choose the following conventions 



r l- X (t)(l-t- 2 )] 2 . (11) 



i' = 1-4^, d2) 

xW = ^lMI f (13) 



G F M% _ 



* (i) = -£*£t=m- (14) 



Among the quantities a, Gf, s^, Mz there are only three independent, and Y z is predicted by the theory as 
well. The phrasing effective Born cross-section means here that we use, besides a (introduced in (0), the 
following input variables: 

4 = 0.23, (15) 

M z = 91.188 GcV, (16) 

T z = 2.495 GeV, (17) 

G F = 1.16637 x 1(T 5 GeV -2 . (18) 

The values are, in a strict sense, related in the Standard Model, and may be determined e.g. by using the 
package ZFITTER JUIH. Here, we took them from iTTTsIl . 

We may now estimate the relevance of the Z-boson exchange to Bhabha scattering in different kinematic 
regions of interest. It is minor at smallest energies where s, |i| << M§, because there xi x ) ~ X /M z << 
l,x = s,t. The strength of the Z exchange amplitude, relative to the photon exchange, becomes at large s,\t\ 
asymptotically: 



G F M§ 



0.3739. (19) 



y/2 87TO 

The other scale of relevance here is the ratio of photon propagators in the s- and t-channels: 

s 2 

1 1 - cos 6< ' 



(20) 



In fact, at meson factory energies, the electroweak Born cross-section agrees with the QED prediction within 
few per mille, and at LEP2 or the ILC within better than 50 %, while at LEP1 or at GigaZ the ratio may become 
bigger than 25; this happens of course only for large scattering angles. At small angles, the corrections may 
safely be normalized to the QED Born cross-section everywhere. The gross features are illustrated in FigureQ] 
for large and small angle Bhabha scattering. For large angles, we show the cross-section ratio separately for 
LEPl/GigaZ and the ILC in Figure [2] We conclude that only for large angles at LEP 1 energies it is better to 
relate the corrections from higher order contributions to the weak Born prediction, while for all other kinematics 
one may use the simple QED Born cross-section. 



III. THE VACUUM POLARIZATION 

Higher-order fermionic corrections to the Bhabha-scattering cross section can be obtained inserting the renor- 
malized irreducible photon vacuum-polarization function, n, in the appropriate virtual-photon propagator, 



q z + io q A + io (f + io 



iX iX ^ q 9 q q ) n ( g ) „2 I 77 - ( 21 ) 
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FIG. 1: Ratio of electroweak to QED Bhabha scattering cross-section at large angles (up) and small angles (down) as a 
function of \fs. 



Here q is the momentum carried by the virt ual p hoton, S — > 0+. The vacuum polarization II can be represented 
by the once-subtracted dispersion integral llll9ll : 



a 2 C°° 
U(q 2 ) =-— dz 

t Jam 2 



ImII(z) 



1 



z q 2 — z + i5 : 



(22) 



where the appropriate production threshold for the intermediate state in II is located at q 2 — 4M 2 . We leave as 
understood the subtraction at q 2 = for the renormalized photon self-energy. 

Contributions to II arising from leptons and the top quark can be computed directly in perturbation theory, 
setting M = rrif in Eq. ( f22l . where mj is the mass of the fermion appearing in the loop, and inserting the 
imaginary part of the analytic result for II. 

We have at one-loop accuracy: 



Im 11,(2) 



(?) 



F, 



m 

m 



f I QJC f e(z-4m})l{^-[3-(3J(z) 



+ epf(z)[3+-L f3f {z) 



4 #, ^ ^ t 

■ 0f(z) - -±r- Lt 



0f (z)]} + O{a 2 ), 



(23) 



3' /v ' 2 

where Qf is the electric charge, Qf = —1 for leptons, Qf = 2/3 for up-type quarks and Qf = —1/3 for 
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FIG. 2: Ratio of ' electroweak to QED Bhabha scattering cross-section at large angles in the energy ranges of LEPl/GigaZ 
(up) and ILC (down). 



down-type quarks, and Cf is the color factor, Cf = 1 for leptons and C/ = 3 for quarks. In addition, we have 
introduced the 9 function, 9{x) — 1 for x > and 9(x) — for x < 0, and the threshold factor, 



L Pf (z) = 
The overall regularization-dependent factor reads as 



in-. 



4(3 2 Jz) 



F r 



ml ir e 



IE 



(24) 
(25) 

(26) 



where fi is the 't Hooft mass unit and -fE is the Euler-Mascheroni constant. 

The inclusion of the 0(e) terms in Eq. (f23T > deserves a comment. These terms might play a role when 
combining Im 11/ with a pole term of another one-loop insertion in a reducible two-loop Feynman diagram. 
The Bhabha-scattering cross section we are going to consider is an infrared-finite quantity, provided one takes 
into account the real emission of soft photons. Therefore, when summing up all contributions, the result does 
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FIG. 3: The one-loop topologies for Bhabha scattering. The gray circle in (c) denotes the vacuum polarization under 
consideration, which may be understood to include fermionic and hadronic one- and two-loop irreducible self-energy 
corrections. 



not show any pole in the e plane and all radiative corrections, including the one-loop photon self-energy, can be 
eva luate d at (9(e). However, we retain the higher e order in Eq. d23l ) for comparing partial results with those 

of eh. 

In contrast to leptons and the top quark, light-quark contributions get modified by low-energy strong- 
interaction effects, which cannot be co mpute d using perturbative QCD. However, these contributions can 
be ev aluat ed using the optical theorem Il20ll . After relating Imllhad to the hadronic cross-section ratio 
i?had CH], 



Imn had (z) 

-Rhad(^) 



3 



Rhad(z), 



hadrons} ; 



(47ra 2 )/(3z) 



(27) 
(28) 



Im n^ad can be obtained from the experimental data for i?had m the low-energy region and around hadronic 
resonances, and the perturbative-QCD prediction in the remaining regions. The lower integration boundary is 
given by M = m n , where m^ is the pion mass. For self-energy corrections to Bhabha scattering at one-loop 
order this was first employ ed in IU411 . Two-loop applications, similar to our study, are the evalua tion of the 
hadronic vertex correction II 1 2 ill and of two-loop hadronic corrections to the lifetime of the muon Il22ll . The 
latter study faces quite similar technical problems to those met here, like the infrared divergency of single 
contributions and the existence of several scales. 

For the fermionic and hadronic corrections to Bhabha scattering at one-loop accuracy, there is only the 
self-energy diagram shown in Fig.[3jc). The two-loop irreducible self-energy contributions have the topology 
shown in Fig.[3jc). One has additionally the four classes of two-loop diagrams shown in Fig.|4]The reducible 
self-energy (Figure H|a)) and vertex (Figure EJb)) topologies are much easier to evaluate than the irreducible 
vertex (Figure 0c)) and box (Figure Hd)) topologies. In fact, only the two-loop boxes were unknown until 
quite recently. 

The two-loop corrections have to be added with the loop-by-loop contributions (the interferences of the 
topologies of Fig. [3]) and with the soft photon corrections. All these terms will be discussed in the following 
sections. 

To summarize this section, the hadronic and heavy-fermion corrections to the Bhabha-scattering cross section 
can be obtained by replacing appropriately the photon propagator by a massive propagator, whose effective 
mass z is subsequently integrated over. Inserting d22b and d27b in (f2T~b we get: 



9^ 



dz 



R(z) 



1 



9^ 



%9v 



q 2 + i 5 3tt J4M 2 z q 2 — z + iS \ q 2 + i 5 

In the following, we will call the massive propagator function in ( |29l > the self-energy kernel function: 



K S E{q 2 ;z) 



1 



7 2 — z + i S 



The weight function R(z) is given by the sum of the non-perturbative light-quark component of Eq. 
the perturbative result of Eq. d23l ). valid for leptons, / = e, /i, r, and the top quark, / = t: 



(29) 

(30) 
and 



R(*) = <adM 



(5) , 3 



E lmU A z ) 



f=e,n,T,t 





(6) 




(<■/) 



FIG. 4: Two-loop topologies for Bhabha scattering with vacuum polarization insertions: reducible self-energy (a) and 
vertex (b) corrections as well as irreducible vertex (c) and box (d) corrections; for the irreducible self-energy corrections 
see Fig. \3\c). 



,(5) 



*») + 



R f (z;m f ) = Q)C f 



J2 Rf(z;m f ), 

/=e,(U,-r,« 




(31) 



(32) 



Compared to ( 123) . we omit here the terms of order O(e). The function i?{ d (z) will be discussed in Appendix 

m 

Corrections related to electron insertions (/ = e) will be discussed separately. For pure self-energy insertions 
(see Appendix lAl. we may consider the electron mass as being small and neg lect terms of order 0(rn\jx), 
x = s, \t\ , \u\. At the expense of that, even the three-loop corrections are known [ 123tl . For two-loop irreducible 
vertex and box corrections, we may either consider m e being finite and treat a two-scale problem (s/m 2 ,, t/ml), 
or we may assume also here m\ « s, |i|, |m|. Instead, for the diagrams with self-energy insertions of other 
fermions /, we will assume rn%, << mi, s, \t\, \u\, but we will make no additional assumption on mi. 



IV. PURE SELF-ENERGY CORRECTIONS 

The pure vacuum polarization contributions to Bhabha scattering form a gauge invariant subset of diagrams. 
So, their numerics may be discussed separately. They can be readily obtained from the tree-level result (O by 
introducing appropriately a running fine-structure constant a(x), where x = s,t, 



da arl 



d9. 



1 

2^ 



>yi{s,t) 



\ a \ s )\ 5 \~2a(t)Kea(s) \- a (t) 



v 2 (s,t) 



St 



and where the running of a is defined as 



Vi(t,s) 



t 2 



0{ml), 



(33) 



a(x) 



1 - Aa(x) ' 



„( 5 ) 



(34) 
HI (see Refs. 1121 



Here Act is given by the sum of the non-perturbative light-quark contribution Aoi ha( 

Il26l Il27ll and references therein for recent developments), a perturbative electron-loop component evaluated 

in the small electron-mass limit, II e , and a fermion-loop term computed exactly, 11/, with / = /i, r, t, 



Aa(x) = Aa| 1 5 a ) d (x)+n e ( a; )+ £ U f (x), 



(35) 



A^Or) 



a x 

7T 3 



dz ■ 



4^) 



K SE (x;z), 



(36) 



with the self-energy kernel function Kse(x; z) < f30b . 

For x < 4m£, Eq. d36] l is well defined. For x > 4ra^, the real and imaginary parts are after a subtraction: 



Re 
1m 



MM 



Aa^^] 



a a; 
^3 



4m 

(5) 
had 



d^ 



(*). 



? (5) 



(5) 



cw-cw 



z (x — z) 



l^dMlo? 



4mJ 



-1 



(37) 
(38) 



The Im Aa^ d (x) coincides with Eq. ( 1271 ). Expressions for the perturbative contributions to the photon 

vacuum-polarization function, 11/ and II e , are available in QED exactly up to two loops [128] and in the 
small electron-mass limit up to three loops Il23ll . For convenience, their explicit expressions are collected in 
Appendix|A] For our analysis, we use the exact results of Eqs. dA3b and (IA4b for fermion loops (/ ^ e), and 
the high-energy expressions of Eqs. ( IA8b , JA9t and ( lAlOb for electron loops. 

In Tables U and |ll] we show numerical values for the various components of Aa of Eq. (f35T > for space- 
like and time-like values of x (t- and s-channel). Note that Aa develops an imaginary part in the s-channel 
above the two -part icle production threshold (see TablelJi. Besides the Fortra n package hadr 5 . f for hadronic 
contributions Il29ll . we employed the Mathematica package HPL II 1 3 OL 1 1 3 1 II and, as a cross check, our Fortran 
routines (see AppendiceslAlandlFii. 



y/s [GeV] 


1 


10 


Mz 


500 


1 loop e 


104.462-24.3245 i 


140.119-24.3245 i 


174.347-24.3245 i 


200.698-24.3245 i 


M 


21.352-24.3060 i 


57.551-24.3245 i 


91.784-24.3245 i 


118.136-24.3245 i 


T 


-0.508 


12.194-24.1724 i 


48.060-24.3245 i 


74.429-24.3245 i 


t 


< 10~ 3 


- 0.007 


- 0.595 


-5.180-29.0633 i 


2 loops e 


0.258-0.0424 i 


0.320-0.0424 i 


0.380-0.0424 i 


0.426-0.0424 i 


M 


0.123-0.0487 i 


0.177-0.0424 i 


0.236-0.0424 i 


0.282-0.0424 i 


T 


- 0.005 


0.118-0.0626 i 


0.160-0.0426 i 


0.206-0.0424 i 


t 


< 10~ 3 


< 10~ 3 


- 0.002 


0.061-0.0876 i 


3 loops e 


0.001-0.0005 i 


0.002-0.0006 i 


0.003-0.0008 i 


0.004-0.0009 i 


hadrons 


-74.420-37.9089 i 


138.850-97.4106 i 


276.213-97.2980 i 


370.744-97.2980 i 


SUM 


51.263- 86.6310 i 


349.324 - 170.3800 i 


590.586 - 170.3997 i 


759.806-199.5505i 



TABLE I: Contributions to Aa in units of 10 4 in th e s-channel (see Eq. Q5t). The real part of the hadronic contributi ons i s 
obtained with help of the subroutine hadr 5 . f Il29tl , the imaginary part follows from the Burkhardt parametrization Il32ll . 




FIG. 5: Hadronic and fermionic irreducible vertex diagrams. The gray circles mark die corresponding one-loop insertions. 



0[°] | xAi[GeV] 


6> = 20 j 1 


6 = 20 | 10 


6 = 3 | Mz 


6 = 3 500 


1 loop e 


77.3512 


113.008 


117.935 


144.286 


M 


3.3069 


30.614 


35.463 


61.727 


r 


0.0148 


1.346 


2.365 


18.804 


t 


< 10~ 4 


< 10" 3 


< 10~ 3 


0.012 


2 loops e 


0.2109 


0.273 


0.282 


0.327 


A* 


0.0260 


0.126 


0.136 


0.184 


r 


0.0001 


0.011 


0.019 


0.097 


£ 


< 10~ 4 


< 10^ 3 


< 10~ 3 


< 10" 3 


3 loops e 


0.0006 


0.001 


0.001 


0.002 


hadrons 


2.6072 


57.830 


71.643 


162.280 


SUM 


83.5177 


203.209 


227.844 


387.719 


6 = 90° j ^s [GeV] 


1 


10 


Mz 


500 


1 loop e 


99.0951 


134.752 


168.980 


195.331 


M 


17.4725 


52.200 


86.418 


112.769 


T 


0.2412 


10.841 


42.746 


69.064 


t 


<10~ 4 


0.003 


0.284 


6.208 


2 loops e 


0.2487 


0.311 


0.370 


0.416 


M 


0.0924 


0.167 


0.227 


0.273 


r 


0.0021 


0.068 


0.150 


0.196 


i 


< 10~ 4 


< 10~ 3 


0.001 


0.021 


3 loops e 


0.0009 


0.002 


0.003 


0.003 


hadrons 


25.0834 


127.219 


256.279 


362.375 


SUM 


142.2363 


492.396 


555.458 


746.656 



TABLE II: Contributions to Aa in units of 10 4 in the t-channel for three values of the scattering angle, 
and 8 = 90°, t = — s sin 2 (8/2). See the caption of Tab.|I]for further details. 



3°, 



20° 



V. IRREDUCIBLE VERTEX CORRECTIONS 



Hadronic and heavy-fermion irreducible vertex corrections are obtained through the interference of the di- 
agrams of Figure [5] with the tree-level amplitude. The contributions from the irreducible vertices are gauge 
invariant by themselves. Their contribution to the 0(a 2 ) differential cross section is given by 



da v 



<in 



(- 

V TV 



a 
2s 



^M ReF2(8)+ !lM F2 ( t) x^M 



t 2 



St 



ReV 2 {s) + V 2 (t)] }+0(m2).(39) 



Here V% summarizes all two-loop fermionic corrections to the QED Dirac form factor, whose computation can 
be traced back to the seminal work of Refs. 113311 and 1 134]. The full result can be organized as 



V 2 {x) = V 2 e{x)+V 2les t{x), 



(40) 



where V 2e denotes the electron-loop component. Closed analytical expressions in the case of electron loops 
at fin ite ra e can be found in Ref. 19511 . In the high-energy limit, compact expressions are available thanks to 
Ref. fl35ll: 



ir ,, 1,3/ ml\ 19 , 2 / m 2 \ 1 /265 ,\, / mt\ 1 /383 , \ , nl 2 , „ n 

After a combination with soft real electron pair emission contributions ( 1DU , the leading logarithmic contribu- 
tions ln (s/m 2 .) get cancelled in i 



10 



Heavy-fermion and hadronic contributions, instead, can be evaluated as in Ref. 1112111 through the dispersion 
integral 



V2rcst(a;) 



dz Ky(x + id; z), 

AM 2 z 



(42) 



where R is given in Eq. OTb and the two-loop irreducible vertex kernel function Ky, in the limit of a vanishing 
electron mass, reads as 



Kx 



{x > z) = \{-\-T x 



3 z \ , ( x 

- H In 

A 1x) V z 



1 / z \ z 



<2 



Li 2 1 



!)]}■ 



(43) 



Here Lia(a;) is the usual dilogarithm and (2 = Li2(l) = 7r 2 /6. The kernel is at the upper integration boundary 
of the order 0(1/ z), the integrand of order 0(1/ z 2 ) so that the dispersion integral is finite there: 



Ky(x;z) 



If 11 

3l36 U 

+ f- 4? 



1 , / N I 13 1,, , 

6 UlnHi)+ 1^288 + 24 ln( ^' 



1 



ln(— u) 



)„•} 



V3600 
At the lower integration bound, the integrand becomes for small z/x: 



for u 



x 
z 



(44) 



Ky(x- z) « 1 { -I - C(2) + J ln(-«) - J ln 2 (- 



l + 2C(2) + -ln 2 (- M ) 



l) 



for u = — 



(45) 



This asymptotic behavior yields at most terms of the order of In (x/M 2 ) if M 2 « x. 

An interesting question is the identification of mass logarithms in case of fermion insertions. Let us rewrite: 



V 2rcst (x) = V£l d (x)+ ]T Q 2 f C f V 2f (a 



(46) 



r(5) 



where V^.^ denotes the non-perturbative light-quark term and V2/ the perturbative contribution of a fermion 
of flavor / 7^ e. Potentially large logarithms arise from parts of the integrand for the z integration which are 
singular at the lower integration bound, z — ► AM 2 , when allowing thereby M 2 to become small. For fermions, 
one has to analyze Rf(z)K\/(x; z)/z in that limit. 

The corresponding analytical integrations may be performed easily after applying the transformation 



Ami 

7 = — 

1 — u A 

thereby getting rid of the square root function in Rf(z): 

R f (z) = C/Q^(3- 
After that transformation, the dispersion integral becomes: 

1 1 



V 2f (x) 



dn 



1 



1 



Ky x + i5: 



Am 2 



(47) 



(48) 



(49) 



From the vertex kernel function Ky (x; z), we have additionally dependences on Ln(— x/z) and on Li 2 (l+x/z). 
Although after the variable change d47b the arguments of logarithm and dilogarithm become non-linear, all the 
integrals may be taken trivially, and we will not go into further details. The result contains L13 and powers of 
logarithms \n n (x/m 2 f) with n < 3. In fact, one will redi scove r in the kinematically interesting ultra-relativistic 
case the formula known from Il35ll and e.g. also from 111 ill : 
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FIG. 6: Hadronic and fermionic reducible vertex diagrams. The gray circles mark the corresponding one-loop insertions. 



1 /3355 19 , n ,\ m , 2 , 

1 C2 — 2 C3 I + 0(m 2 f ). 

6 V 216 6 S V V /; 



(50) 



The same soft- real pair cancellation mechanism as described for electrons works also for heavy fermions, and 
the leading logarithmic powers In (s/rrA) will get cancelled in the cross-section. This is of physical relevance 
if the soft pair emissions remain unobserved. In our numerical studies, we will, conventionally, include the soft 
electron pair emission cross-section, but not that fo r hea vy fermions or hadrons. For further det ails s ee Section 
iDl and some numerical results were presented in jl36ll . where we used the parameterization Il32ll with flag 
setting I PAR = 1. 

We just mention that the transformation J47] ), when applied to the simple one-loop self-energy kernel (|30l l. 



K S e(x;z) 



1 



4m 2 f /x 



Arrif/x 



(51) 



gives a rational integrand for the it-integration, and one gets as a result a function at most linear in \xi(s/mV). 
For the explicit expressions see Equations (IA3b (constant term in e) and dA8b . 



VI. INFRARED-DIVERGENT CORRECTIONS 



There are various origins of heavy-fermion or hadronic infrared divergent cross-section contributions of 
order 0(a 4 ): 

• Factorisable diagrams with one-loop vertex or box insertions 

• Irreducible two-loop box diagrams 

• soft real photon corrections 

The sum of these corrections is gauge-invariant and infrared finite. 
We will consider five classes of contributions: 

(a) Interference of Born diagrams with reducible [vertex+self-energy] corrections of Fig.[6j 

(b) Interference of one-loop vertex and self-energy diagrams, both of Fig. [3j 

(c) Interference of one-loop box and self-energy diagrams, both of Fig. |3j 

(d) Interference of real soft photon emission diagrams, one of them with a self-energy insertion; 

(box) Interference of Born diagrams with two-loop box diagrams of Figure|7] 

For ease of notation, in the following we collect the overall dependence on a and rewrite the factorizing 
contributions of class i,i = a,...,d: 



da 



fact. 



dQ 



2 a 2 da] 



fact. 



dn ' 



(52) 



and analogously for the two-loop boxes. In addition, we define 

s 



)))- 
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(53) 
(54) 



and introduce short-hand notations for those kinematic factors which appear more than once in the following 
formulas: 



(_ .-■/■■ 
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s 2 st r 
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C T | ■ 
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(l + (l-r) 2 ) 



(55) 



A. Factorisable corrections with vertex or box insertions 

The infrared-divergent fa ctoris able heavy fermion and hadronic corrections for m 2 << M 2 , s, \t\, \u\ can be 
readily obtained from Ref. ill 1 ill by replacing the photon vacuum-polarization function in the s- or ^-channel 
with the dispersion integral 



o 



n(x) = Aa(x) = - I(x), 



I(x) 



dz R(z) 



3 Jam 2 z x — z + id' 



where Aa(x) is given in ( f35l > and R in QTI ). 

We begin with the reducible vertex corrections (a). From Eq. (3.8) of Ref. ill 1 ill we derive: 



rf(T fact. 
dfl 



Fe 
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(Ar (l-ln(sj)ReI{s)- Trim I(s) + B r ln(s) +ln(r) - 1 I(t)\ 
+ i | A. In 2 (J) - 8 Ca - (A - 2 C* r ) ln(l) + 2 ^ r - C r \ } Re /(s) 
2A r ln(l) - A + 2 CV 7rlm/(s) - - |s r ln 2 (J) + ln 2 (r) 
E r -2B r ln(r) ln(s) - S r ln(r) - 2 S r < 2 + 8F r | J(t), 



1 r 

2 



where the normalization factor i 7 ^ is given in Eq. ( I26l l. It appears here in the combination 



F, 



1 , 
In 



1h(tt) - 7E + 0(e). 



(56) 
(57) 



(58) 



(59) 



In strict analogy, the interference of the one-loop vertex diagrams of Fig ure [3] (a), with the vacuum-polarization 
diagrams of Figure[3](c) can be extracted from Eq. (3.26) of Ref. Ill ill : 



dn 
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[ A (l - ln(J)) - D r ln(r) ReJ(s) - C r 7rlm/(s) 
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A (ln(J) - l) + G r ln(r)] I(t) } + ^ { A r ln 2 (s) 



-2 



1 - 4 r A - 4 r" 
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}7rIm/( S )-i{Aln 2 (.s) 



£ r - 3 (l + rA-) ln(r)l ln(s) + — (l + r A-) ln 2 (r) - — x 



r 
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(60) 



Final ly, th e contributions from the one-loop box diagrams of Figure [3] (b) may be derived from Eq. (3.28) of 
Ref. mill: 



da] 



fact. 
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+ -(2D r + r) ln(l-r) + ^ (l - r) ln 2 (r) + - (l-2r) ln 2 (l - r) 



1 

27 
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x ln(r) + - 3(l-2r) +4r 2 ln(l - r) + -^ (r A + 1 + 2r 2 ) } Trlml(s) 

+ { [ A ln(r) + A ln(l - r) ] ln(s) + — ((J r + 2) ln(r) - 



1 
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r A 



2 (1-r 

1 
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ln(l-r) + — (5-4r] In 2 (r) + — (2 - r) ln 2 (l - r) 



2rA + 21-r 
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2r 



•)c a }/(*). 



(61) 



All three types of corrections are infrared divergent. The vertex diagrams contribute leading electron mass 
singularities of the order In (s/m 2 .), while for the factorisable box diagrams the leading order is ln(s/m 2 ). In 
addition, the self-energy insertions I(x) yield a dependence on In^/mi), in case mi is small compared to s. 



i 2 f), in case ml 
This may be most easily seen from the e-independent terms in (|A3t . So, we collect here at most terms of the 

4) 



order In 2 (s/m 2 ) ln(s/m 2 f ). 



B. Soft real photon emission 

In order to obtain an infrared-finite quantity, we take into account the interferences of diagrams with real 
emission of soft photons from the external legs, where one of the diagram has a vacuum-polarization insertion. 
The anatomy of these real corrections is exemplified in AppendixICl where the soft photon factor is shown both 
for non-vanishing elec tron m ass m e and in the ultra-relativistic approximation. The result may be also read off 
from Eq. (4.4) of Ref. Ill ill and reads as 



da'. 



fact. 



da. 



4,1 

fact. 



</o 



da\ 



2w 
dtt ' [Jsj dfl 



d,1 

fact. 



(62) 



where u is the maximum energy carried by a soft photon in the final state. We obtain 
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fact. 
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ln(s) + ln(r) - ln(l - r) - 1 A r Re I(s) - B r I(t) 
2 J - ln 2 (s) + ln(s) fln(r) - ln(l - rj\ + - ln 2 (r) - - ln 2 (l - r) 
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FIG. 7: Irreducible box diagrams. The gray circle denotes the hadronic or fermionic insertions. 



ln(r) ln(l - r) - 2Li 2 (r) - (2 A r ReI(s) - B r I(t) 



D r 



ln(s) + ln(r) - ln(l - r) - 1 



ReiYs 



)+'(*)]} 



da 



fact. 



(in 



ln(s) + ln(r) - ln(l - r) - 1 A r Re I(s) - S r /(*) 



(63) 



(64) 



Again, the infra-red divergency is contained in the factor F e /e, and the mass singularities are at most of the 
orders hx(x/rrA), x = s,t, and In (x/rn^) for the ^-independent part and hi{x/m1) for the w-dependent part. 

C. Two-loop irreducible box corrections 



From the technical point of view, the two-loop irreducible box corrections of this section, represented by the 
three box kernel functions, are the main result of the article. Their contributions to the Bhabha-scattering cross 
section arise from the interference of the diagrams of FigureQwith the tree-level amplitude and can be written 

as 



d(Thoi 



da] 



box 



2 a 2 
— 2 , 

4s V 



dfl 
(ReA s 



Re A t 

t 



(65) 



Here the functions A s and A t contain the interferences of box diagrams with the s-channel and t-channel 
tree-level diagrams and can be represented through three independent form factors, evaluated with different 
kinematic arguments: 



A s = B A (s,t)+B B (t,s)+B c (u,t)-B B (u,s), 
A t = B B (s,t) + B A (t,s) - B B (u,t) + B c (u,s). 



(66) 
(67) 

In addition, note that in Eq. ( l65l l we have collected an overall factor 1/4, coming from the sum over the spins, 
and a factor 2, taking into account the fact that the contributions generated by the diagrams (la), (2a), (3a) 
and (4a) are equivalent to those of diagrams (16), (26), (36) and (46) of Figure|7] Finally, the correspondence 
among the form factors of Eq. d66l i and the diagrams of Figure|7]reads as follows: 

B A {s,t), diag. 1 x trcc t => B B (s,t), 

B B (t,s), diag. 2 x tree t =$• B A (t,s), 

Bc(u,t), diag. 3 x treet => —B B (u,t), 

—B B (u,s), diag. 4 x trcc t => B c {u,s). (68) 



diag. 1 x trec s 

diag. 2 x tree s 

diag. 3 x tree s 

diag. 4 x tree s 
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FIG. 8: The one-loop master integrals with an additional mass scale M — y/z for the dispersive two-loop box evaluation. 



We evaluate the three form factors Bi using dispersion relations and computing thereby the convolution of 
the hadronic or fermionic cross-section ratio R with three kernel functions A',, 



, , f°° R(z) , 

Bt{x,y)=\ dz Ki{x,y;z), 

Jam 2 z 



(69) 



where R has been introduced in Eq. (|3TT i. and the kernel function are to be calculated. For positive x or y, one 
has to replace x — > x + iS or y — > y + iS. 

The self-energy insertion is represented by a dispersion relation, thus replacing the one-loop photon propa- 
gator by a massive effective propagator as in Eq. ( 1291 . This procedure reduces the evaluation of the two-loop 
diagrams to one-loop complexity with a subsequent dispersion integration. Employing standard techniques, 
together with the Mathematica packages AMBRE 11373 and MB II138I1 . for a reduction of one-loop integrals to 
scalar master integrals, the kernel functions have been finally expressed by eight one-loop master integrals 
MU)( x ,y;z), 



K l (x 1 y;z)=F e Y,C { l 3 \x,y- 1 z)M^(x 1 y;z), 

3=1 



(70) 



where F t is the usual normalization factor of Eq. (1261 1. and C\ are rational functions of the kinematic in- 
variants, of the space-time dimension d, and of the two masses m e , z. The master integrals M^) are shown 
in Figure |8] and analytical expressions for them can be found in Appendix iBl Due to their length, we do not 
reproduce here the explicit (exact in m e and d dimensions) right hand side of (l70l . but refer for them to the 
Mathematica file at the webpage llll7ll . 

In the small electron-mass limit we obtain the two-loop box kernel functions: 



z (x + y) + x" 



3(y-z)*. e 

+ 2x(x + ij)+ y 2 

-y(2x + y) 



z[ \-y 

y 



In 



X 



ln 2 f-^ 



z 2 + 2zx 



2z z [- + l)-z[ h6x + 5y 



U 



!) 



+ x (x + Ay) + 'iy 2 lnf-— \ + z 2 - 2 z (— + x + y) 



16 



+ 2x (x + y) + y z 
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x In' 
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x In I — 






(73) 



- 1 Li 2 - - 2 x z Li 2 1 

These kernel functions are reproduced in Mathematica files at the webpage Ill7ll as functions KA, KB, KC 
and KAexp, KBexp, KCexp. 

The two-loop box kernel masters (T7TT > to (l73l are evaluated in the Feynman gauge; they are infrared divergent 
and contain collinear singularities in m e . 

After inserting Eq. ( fTTI i. Eq. d72l and Eq. ( f73l in Eq. ( |69l l, we derive the total contribution to the cross section 
generated by box diagrams. Collecting powers of a, we write 
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where the integrand functions are given by 
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The functions 7i (z) to -Z3 (z) are reproduced as functions Ii,l2, I3 in a Mathematica file at the webpage Ill7ll . 
Note that, after assembling all irreducible box diagrams, their total contribution is free of collinear divergen- 
cies in m e because ln(mg) vanishes in the combination 
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ln(7r) - ln(^ 



0(e). 



(78) 



This fact might be observed already for any sum of single pairs o f dire ct and their related crossed box diagrams, 
which is gauge-independent and free of collinear singularities Il39ll : from d68l ) and Figure [7] one selects e.g. 
the following ones: 



K B (t,s;z) - K B (u,s;z), 
K A (s,t;z) + K c {u,t;z). 



(79) 



In the limit nil « s,\t\,\u\, the z-integration over the Ii(z),i = 1,2, develops mass singularities from the 
lower integration bound: 



dz K SE (y;z) A(x,y, z) + B(x,y) In 



(80) 



4M 2 



where A, B are regular for z — > 0. It follows immediately that the irreducible box diagrams yield terms of 
the order of at most In (s/m 2 f), because A joins, after integration, terms with a behavior like a one-loop self- 
energy, and B joins terms with one order more in the logarithmic structure. This has been discussed already in 

CM. 

The residual infrared-singular part of the box cross-section is: 
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The function I(t) (see Eq. ( |57] >) stems from diagrams with a vacuum polarization insertion in the t-channel, 
and I(s) from insertions in the s-channel. One may wonder which of the other infrared divergent parts are 
needed to compensate the double-box divergency (in the gauge chosen here). This may be exemplified by 
collecting all the IR-divergencies of the diagrams with a vacuum polarization insertion I(t) in the t-channel; 
for the others, quite analogue arguments hold. From Sections TVI Al and IVIBI we may extract such terms. There 
are the following divergencies due to vertex diagrams: 
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The reducible box diagrams are (in the curly brackets) free of electron mass singularities, also in the terms not 
shown here. They depend also on u: 
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For the soft real terms, we refer to Appendix [C] and may distinguish between initial and final state corrections 
(which are equal) and the initial-final state interference: 
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It is now easy to see that the IR-divergency of the double box diagrams, , being proportional to ln(— u/s), gets 
completely cancelled by the sum of the reducible box diagrams and the interference part of soft bremsstrahlung. 
Although, the latter introduce to the sum an IR-divergency with ln(— t/s), and this gets cancelled the reducible 
vertex diagrams, thus introducing an IR-divergency with ln(s/rrig), which will be cancelled finally by the initial 
and final state soft corrections. The lesson is: a sensible, infrared safe cross-section contains the complete sum 
of all the single IR-divergent diagrams, or no one of them. 

Despite of that, an isolated treatment of the pure self energies or of the irreducible vertex corrections is 
possible. 

Finally, we just mention that the analytical integrations over z may be performed following the hints in 
Section |V] 



D. Kernel functions for the infrared safe sum 



We are now in a state to evaluate the net cross-section contribution from the various infrared divergent terms 
of Sections lVI Al and lVI CI We have seen that they have to be treated together. The sum of the box contributions 
of Eq. ( F74l > with all infrared-divergent factorisable corrections, given in Eq. d58l , Eq. d60l l. Eq. d6lT l and Eq. 
is infrared-finite and can be cast in the following form: 
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The lower bound is AM 2 = 4m 2 for hadrons and AM 2 = 4m 2 f for fermions /. The auxiliary functions Fi(z) 
are given by 
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The 13(2) is defined in d77b . For < s < 4 M 2 we can write 
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For s > 4 M 2 , we have to perform some subtractions in order to make the formulas explicitly stable around 
z = s, and at the time retain the sufficiently fast vanishing of the integrand at z — > 00: 
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In the limit ml << s, |i|, |u|, the z-integration over the Fi(z), % = 1,2, develops mass singularities from the 
lower integration bound: 



dz K SE (y;z) 
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(94) 
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where A, B, C are regular for 2 — ► 0. It follows immediately that the sum of all infrared divergent diagrams 
yield terms of the order of at most In (s/ml) and ln(s/m 2 ) ln(s/m?), because A joins, after integration, 
terms with a behavior like a one-loop self-energy, B joins terms with one order more in \n{s/mV) and C goes 
to gethe r with at most ln(s/m 2 ) h\(s / ml); there are no cubic logarithms here. This has been discussed already 

mini. 

Further, for the numerical evaluation, the functions F\, F 2 and F 3 are replaced for z — > 00 by their asymp- 
totic values: 
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VII. NUMERICAL RESULTS AT MESON FACTORIES, LEP/GIGAZ, ILC 



We begin with numerical results for Eq. d87b . multiplied by the overall factor (a/ir) 2 a 2 / s. The expressions 
contain the contribution of irreducible two-loop boxes, summed up with reducible two-loop vertex and loop- 
by-loop diagrams, and combined with soft-photon emission. They are called here 'rest' from electrons, muons, 
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6[°] | Vs [GeV] 


6» = 20 1 


e = 20 | io 


6 = 3 Mz 


9 = 3 | 500 


QED Born 


214.903 


2.14903 


53.0348 


1.76398 


weak Born 


214.903 


2.14930 


53.0376 


1.76390 


QED Born, running 


218.559 


2.23814 


55.5353 


1.90910 


vertices [/i+r+hadr.] 


-0.001086 


-0.00022513 


-0.007982 


-0.00129296 


vertices [e] 


-0.102787 


-0.00325449 


-0.092546 


-0.00574577 


soft pairs e + e~ 


0.130264 


0.00403772 


0.112763 


0.00685890 


rest: e 


0.235562 


0.00497834 


0.135650 


0.00672652 


M 


0.009518 
-0.017214 


0.00135040 
0.00134282 


0.040792 
0.040688 


0.00287809 
0.00287795 


r 


0.000074 

X 


0.00005385 

X 


0.002706 
-0.009610 


0.00087639 
0.00083969 


hadr. 


0.008642 


0.00269490 


0.087618 


0.00810781 



TABLE III: Numerical values for the differential cross section in nanobarns at scattering angles 9 = 20° and 6 = 3°, in 
units of 10 2 . Concerning the finite remainder, containing irreducible box diagrams, we show for each fermion flavor the 
result obtained through the dispersion-based approach (first line) and the one coming from the analytical expansion (second 
line), neglecting 0(m't/x), where x — s, \t\, \u\. When m? > x, the entry is suppressed. 



tau-leptons, and from hadrons. The top influence was also considered but comes out so marginal that we don't 
discuss it. The results are summarized in Table [III] and Table [IV] for small- and large-angle scattering and a 
variety of energy scales. We do not discuss the isolated irreducible two-loop boxes because this would become 
more convention-dependent. Note further that in these tables the dependence on the maximal energy of the 
soft photons is switched off by setting u) = y/s/2 (an analogous consideration holds for the soft pairs e + e~). 
For comparison, the tables also contain entries with pure QED Born, QED Born with running coupling, and 
effective weak Born cross-sections, as well as contributions from: electron vertex insertions and soft e + e~ 
pairs (with a quite small sum of them); the su m of heavy fermion irreducible vertices. The hadronic results 
have been obtained using the parametrization lll32ll with flag setting I PAR = and implementing narrow 
resonances as described in Appendix[E] 

We see that the two-loop corrections from electron insertions (the so-called Nf = 1 corrections) are the 
largest, and the second-largest ones are the hadr onic corrections. The tables also demonstrate that the approxi- 
mation m 2 f « s 



f 



\t\, \u\ as applied in e.g. Ill ill works well in the regions where this is expected. 



nA [GeV] 


1 


10 


M z 


500 


QED Born 


466537 


4665.37 


56.1067 


1.86615 


weak Born 


466526 


4654.16 


1238.7500 


0.92890 


QED Born, running 


480106 


4984.83 


62.9027 


2.17957 


vertices [/i+r+hadr.] 


-16.351 


-2.0437 


-0.125208 


-0.0104275 


vertices [e] 


-477.620 


-12.3010 


-0.298589 


-0.0155751 


soft pairs e + e~ 


648.275 


16.0690 


0.376531 


0.0191990 


rest: e 


807.476 


14.5277 


0.270575 


0.0119285 


M 


160.197 
152.890 


6.0819 
6.0809 


0.147046 
0.147046 


0.0072579 
0.0072579 


r 


2.383 

X 


1.3335 
1.0739 


0.075268 
0.075214 


0.0045713 
0.0045712 


hadr. 


232.674 


16.0670 


0.469944 


0.0246035 



90°, in units of 10~ 4 . 
A more detailed picture of the relevance of the fermionic and hadronic two-loop corrections may be got from 



TABLE IV: Numerical values for the differential cross section in nanobarns at a scattering angle 
See the caption of Tablelllllfor further details. 
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figures [9] to [14] where we show the cross-section ratios 

1Q 3 ^NNLO^ (9g) 

where cLoq is the effective weak Born cross-section at ^s = Mz, 500 and 800 GeV, and the QED Born cross- 
section elsewhere. So, the figures show just the relative size of the corrections in per mille. For a comparison, 
we show also the pure photonic corrections. The licNNLO is here the net sum of all the terms discussed arising 
from a fermion flavor (e or /i) or from the hadrons. In case of electrons, we add also the real pair correction. 
The total non-photonic term includes als o the r and top quark contributions. For hadrons, we decided to use 
the parameterization i?i ia d,i as given in lll32ll with parameter IPAR = 1. We applied also numerics with a 
combination i?i ia d,n of several adjusted pieces valid at different scales, as explained in Appendix|E] In Figures 
l9l and [TT1 it is seen that the predictions with i?i ia d,i and i?i ia d,n are quite close to each other. Because we 
did not get a stable numerics over all the parameter space with i?i ia d,n, we decided not to use it for the final 
determination of the physical results until we have a better understanding of its behaviour. 

At a meson factory with t/s ~ 1 GeV (Pig-O the heavy fermion effects are below 0.5 per mille and are thus 
certainly negligible. At y/s i=a 10 GeV (Fig.flOl). electron and hadron corrections amount to 2 to 5 per mille and 
might play some relevance. At the higher energies, we have to consider small angles and large ones separately. 
The hadronic corrections amount to up to 4 per mille at LEPl/GigaZ and 20 per mille at ILC energies at large 
angles, while at small angles they stay well below 5 per mille. For ^fs = 500 GeV this is exemplified in Figure 
[131 and from the tables one may read exact values at 9 = 3 degrees: for the infrared-finite remainder containing 

.had 

box diagrams, at LEP/GigaZ it is J eat =1.65 per mille, and at ^/s = 500 GeV the corresponding value 

becomes 4.6 per mille. Everywhere, the pure photonic corrections are the largest one, followed by the Nf = 1 
corrections. This is, of course, due to the small electron mass producing large logarithmic mass effects and is 
extensively discussed in the literature. 



VIII. SUMMARY 

The NNLO effects of heavy fermions and hadrons on the Bhabha cross-sections are accurately known now 
and the determination of QED two-loop corrections is completed. For each of the corrections there exis t 
several independent calculations. Quite recen tly, a seco nd determination of the h adronic corrections in Il40ll 
fully confirmed our results as presented in IUl lll2l[l36ll and at our webpage 1 14- ill . We ind eed c hecked, when 
preparing this longer write-up of our results, th at, w hen using the same parameterization lll32ll . all the digits 
shown in our Table s agr ee with those shown in lll40ll (see Tabs. Hill and llVb . The nume rical diffe rences which 
were mentioned in Il40ll were due to a different choice of the parameter IPAR in Ill2ll and ll40ll . 

Summarizing the numerical discussion, it is quite obvious that for measurements aiming at an accuracy at 
the per mille level it is crucial to take the heavy fermion and hadron contributions into account. A detailed 
conclusion for a specific experiment evidently depends on the experimental set-ups and will deserve the use of 
a precise Monte-Carlo program. 

Finally, we would like to mention that, in pure QED, not all of the contributions have been determined so 
far. It would b e qui te interesting to know also the influence from the so-called radiative loops. This problem 
was treated in Il42ll . but so far without account of the radiative loop diagram, which include e.g. radiative 
boxes with the need of knowledge of five-point functions. Also here, final conclusion will be made only with a 
precise Monte-Carlo program. 

As a third field of future improvement we like to mention the complete treatment of electroweak two-loop 
corrections to Bhabha scattering. As already said there exists some literature on that subject. The leading 
NNLO weak corrections due to top quarks have been determined long ago in flltl . This was considered as a 
satisfactory approximation for LEP1 and implemented e.g. in the packages ZFITTER 18111 and in th e pro gram 
family KORALZ |66|, KKMC HEUl, BHLUMI JH], BHWIDE J42J]; see also the workshop report Il44ll . An 
improvement of that might become necessary for large angle scattering at the ILC. This might be done similarly 
to the recent implementation of weak two-loop corrections for muon pair production in ZFITTER v.6.42 18111 . 
based on original work described in ILZllOl and references therein. 
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FIG. 9: Two-loop corrections to Bhabha scattering at y^s = 1 GeV, normalized to the QED tree-level cross section. 
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APPENDIX A: ANALYTIC RESULTS FOR THE FERMIONIC VACUUM POLARIZATION 

The contribution of a fermion of flavour / to the irreducible renormalized photon vacuum-polarization func- 
tion n, introduced in Eq. (fJTJ, can be written in pure QED as 



n/(s 2 ) = E ©> 5 Qfc f nf{<f) + o{a% 



(Al) 



where Qf is the electric-charge quantum number, C/ is the color factor and the normalization factor F e is 
defined in Eq. ([26}. 

For our purposes we need both the n = 1 and n = 2 terms up to O(e ). However, since some components 
of the infrared-finite differential cross section show single poles in the e plane, we find useful to consider also 
the 0(e) part of the one-loop photon self-energy for intermediate checks of the results. 

Both expressions can be written in a compact form innoducing the variable 



J-q 2 + 4m 2 f - 

J '-q 2 + Am 2 f + 



(A2) 
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FIG. 10: Two-loop corrections to Bhabha scattering at y/s = 10 GeV, normalized to the QED tree-level cross section. 

The results can be found in Appendix A of Ref. lll45ll and at the webpage lUl7ll . In the space-like region 
— oo < q 2 < 0, it is < x < 1, and one gets a real vacuum polarization: 
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(A4) 



For the time-like region, we have to perform an analytical continuation to q 2 > 4 m 2 ? by setting q 2 —* q 2 + i S 
in Eq. (1A21 i. Now, the conformal variable x develops a small positive imaginary part and it is —1 < Rex < 0. 
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FIG. 11: Two-loop corrections to Bhabha scattering at y/s = Mz, normalized to the effective weak Born cross section. 



In order to derive Im II of Eq. d23l ), we may introduce an auxiliary variable y: 



q 2 — ./q 2 — 4m 2 f 



q 2 - 4 m 2 f 



(A5) 



and observe that x = —y + i8, with y = for q 2 — > oo and y = 1 for q 2 = 4 nii. With these conventions, it 
becomes evident for Eqs. dA31 l and dA4l ) that Li2(±a:), Li3(±a;), and ln(l + x) stay well-defined, and one has 
to take care about ln(a;): 



ln(x) —* ln(— j/ + i (J) = ln(y) + J7r. 



(A6) 



Of course, one may perform the evaluations with complex variables either. 

The contribution of electron loops to the irreducible renormalized photon vacuum-polarization function II 
of Eq. (1211 in the small electron-mass limit is available in pure QED up to three loops, 



ne(9 2 ) = E(-)" n ^(? 2 ) + o(« 4 ). 



n=l 



(A7) 



The one- and two-loop contributions can be obtained by expanding Eqs. dA3t and ( IA41 > and neglecting terms 
suppressed by positive powers of the electron mass. The three-loop component, (we d o no t include double- 
bubble diagrams with two different flavours), can be found in Eqs. (7) and (9) of Ref. lll23ll . The results for 
q 2 < are: 
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FIG. 12: Two-loop corrections to Bhabha scattering at y/s = 500 GeV, normalized to the effective weak Born cross section. 
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The continuation to q 2 > is again obtained by the replacement q 2 — ► q 2 + id. 



APPENDIX B: MASTER INTEGRALS FOR THE BOX KERNEL FUNCTIONS 

The three kernel functions for irreducible box diagrams of Figure [7] may be found at webpage [117] with 
their exact dependences on m e and on e. They are expressed by eight master integrals, which were evaluated in 
the limit m 2 « z,s,\t\,\u\. The master integrals of Eq. dTOl i, for x = s and y = t, are evaluated to the power 
in e needed here: 
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FIG. 14: Two-loop corrections to Bhabha scattering at y/s = 800 GeV, normalized to the effective weak Born cross section. 
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For M W and M^ 2 \ results are needed up to 0(e), since, after the reduction procedure, both coefficients c\ 
and cj , for i — A,B, C, include terms C(e _1 ). For all other basis integrals, O(e ) results suffice. Note that 
for M^ 1 ) (tadpole), M ' 3 ' and M' 5 ) (no dependence on m e , apart from the normalization factor N) results are 
exact. In other cases, the order of the expansion in m e depends on the coefficients c\ . For example, we have 
c\ ' = 0(mj 2 ), and we compute A/' 2 ' up to 0(mP e ) (note the overall factor m 2 in Eq. d65ll). In contrast, we 
have c\ = O(m e ) and we do not need M^ up to O(m e ). 



APPENDIX C: SOFT REAL PHOTON EMISSION 



The leading order contributions to the soft real photon corrections 

e~ (pi) + e + (p 2 ) -> e~ (p 3 ) + e + (p 4 ) + 7 ( fc ) 
to the Bhabha cross section (f2]i are contained in the factor F so f t : 



do 



LO 



do~o a 



F so ft(w, s,i,m 2 ), 



with oj being the upper limit of the energy of the non-observed soft photons: 

E 7 e[0,w]. 



(CI) 



(C2) 



(C3) 
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The u) has to be chosen as small as to guaranty that the emitted photon does not change the kinematics of the 
process (|T). The NLO radiative cross section with 0(a) vacuum polarization insertions is: 






^!(!!l%i) R en«(.s) 
s [ s z 

Vi(t,S 



v 2 (s,t) 



St 



Re 



n (1) ( s ) + n (1) (t) 



ReIlW(M (-J F BO ft(w,*,t,m»). 
The result for the soft photon factor is split into initial and final state radiation and their interference: 



where 



F so{t (u), s,t,rn\) = S ini + 5 int + S Rn , 



S iui = {Ql + Ql)F 11 +Q 1 Q 2 F 12 
= 2-Fii + F12, 

*lnt = (QlQ 3 + Q2Qi)F 13 + (Q1Q4 + Q 2 Q 3 )F 14 

= 2F 13 + 2Fi 4l 

San = (Ql + Ql)F 33 + Q 3 Q A F 3i 

= 2F33 + F34. 



(C4) 

(C5) 

(C6) 

(C7) 
(C8) 



Each of the terms in Eqns. dC6b to ( IC8b exhibits the radiating particles - a factor QiQj marks the emission of 
the photons from particles with momenta pi and p.,; Of course, it is QiQj = 1 here. Since the initial and final 
state particles have equal masses, it is additionally: 



F 33 — i*ii, 

F34 = F\ 2 . 



(C9) 
(CIO) 



So, it will be: 



F soh (uj, s, t, m 2 e ) = 4F n + 2F 12 + 2F 13 + 2F 14 . 



(Cll) 



The evaluation of F so ft follows standard textbook methods (see e.g. for details in Sec. (4.3) of Ill46ln . The 
exact result for the soft radiation functions is, for d = 4 — 2e: 
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We use the abbreviations: 



(C12) 
(C13) 

(C14) 
(C15) 

(C16) 



P = y/l - 4m 2 /s, 



(C17) 



[3 t = Vl - 4m 2 /i, 
T = 2m 2 - t, 
A^ = \/T 2 - Am 4 , 



[3 U = y/l - Am 2 /u, 
U = 2m 2 - u, 



Xu = \JU 2 - Am 4 . 
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(CI 8) 
(C19) 
(C20) 
(C21) 
(C22) 
(C23) 



Our kinematics fulfills here s + t + u = 4m 2 , and it is T, U > 0. If necessary, the logarithms and dilogarithms 
may be analytically continued with the replacement 



s — > s + te, 



e.g. 



Li 2 



2/3 



= — Lio 



-L l2 (l)--lB ^ 



In the limit of small electron mass m e , this simplifies considerably (s = s/m 2 ): 

Fu = A e + lln(s), 

Fia = -2A £ ln(s)-iln(s) 2 -2C 2 , 



F13 = -2A e ln 
F 14 = 2A £ ln 
Finally, the divergent part is: 



u \ 1 



iln(5) 2 -2C 2 -Li 2 (-^ 



/))- 



-ln(,s) 2 + 2C 2 + Li 2 (-- 
2 \ u 



F 




/2w 




- In (I) 


-In — 


e 




\V« 



Taking all the terms together, we obtain: 

"F 



F so ft(w,s,i,m e ) 



! (|) - 2 lr ( ^ ] 
V \fs) 



-21n(s) + 2-21n 



(s) 2 - 4( 2 + 21n (I) + 2Li 2 ( -- J - 2Li 2 {-— 



This expression agrees, of course, with e.g. Eq. (4.5) of ill ill . 



(C24) 
(C25) 

(C26) 
(C27) 
(C28) 

(C29) 
(C30) 



(C31) 



APPENDIX D: REAL FERMION PAIR OR HADRON EMISSION 



The numerical influence of the virtual corrections gets modified by the non-observed emission of real pairs 
of electrons or other fermions, or of hadrons: 



dcr rcal da a 2 

dn ~~ In^ 2 



[S e + S f + 5 



had'] 



(Dl) 



The real pairs or hadrons give non-singular contributions and depend, in the simplest configuration, on an 
energetic cut-off D on the invariant mass of the non-observed pair or hadrons F rea i, and of course also on the 
production threshold 2M. 

There are two basically different situations. In case 4M 2 << s, |t|, \u\, one may additionally choose 2M < 
F r eai < DEbc&m << Fboam (remember Fbeam = y/s/2), and observes a logarithmic dependence of the 
cross-sections on the two parameters M, D. In the other case, assuming M >> m e but otherwise arbitrary, as 
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it is done in the present study if not stated differently, the concept of soft pairs becomes senseless and one has 
to evaluate the pair and hadron emission cross-section numerically with MC methods. 

For completeness and because of the numerical importance, we will include the soft pair emission contri- 
butions for electrons, which is by far the biggest one. For this case, analytical expressions with logarithmic 
accuracy are known from [45]: 
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The parameter D has to fulfill: 



2m e << DE h 



l±cos(9 



«E h 



v = t, u. 



(D2) 

(D3) 
(D4) 
(D5) 
(D6) 

(D7) 



From the sum of dDU and ( |39t , the compensation of the leading mass singularities (contained here in the 
L^,Lf, L\ terms) in the cross-section becomes evident. 



APPENDIX E: THE CROSS-SECTION RATIO i? h ad 



The numerical values of the irreducible two-loop corrections depend crucially on i?had(s) as defined in 
\, while the reducible corrections may be evaluated with one of the publicly available parameterizations of 
II(g 2 ) (see d22b). Unfortunately, we did not find an actual, publicly available code for i?had(s) that covers the 
complete integration region from the thre shold at s = 4M% to infinity. In our short communication \m, we 
used the Fortran rout ine of H. Burkhardt Il32h . This parameterization dates back to 1986 and was used for 
the numerics in ill 2 ill , and it was available by contacting the author Il32ll . The Fortran file is made available 
at our website Ill7ll . I t is t o be expected that current hadronic data would not induce changes compared to 
the parametrization of lll32ll of more than about 10%. This would be tolerable in view of the smallness of 
the irreducible two-loop contributions in our analysis. For the numerically much more sensitive reducible 
contributions, the running coupling a em is needed, and implementations of that are publicly available, e.g. the 
Fortran package hadr5 . f at (129]. 

For the present study, we improved our numerical basis for the evaluation of the irreducible vertex and box 
contributions by combining packages for the evaluation of i?had(s) in different kinematical regions: 

2. 



(El) 
(E2) 



(A) From threshold at s = 4m^ to s = 0.03 GeV^: We follow Section 8.1 of 114 



i?had(s) = i? 7r + 7r -(s) = - I 



Ami 



2\ 3/2 



F n {s) = 1 + 1.879 



( s 



VGeV 



I^WI 2 , 



3.3 



GcV 



-0.7 



GeV 2 



The a bove is based on a fit to e + e 
lll48ll are also taken into account. 



data whose results are shown in Table 3 of 114711 : space-like data 



33 



(B) From s = 0.03 GeV 2 to s = 10000 GeV 2 : Use of subroutine 114. 

(C) Above s = 10000 GeV 2 : Use of subroutine rhad . f v . 1 . 0, published in 115 

In Figure (TT3T > we show the i?had resulting from our Fortran implementation for the regions (A) to (C) as 
described above. 
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FIG. 15: The implementation o/i?had used for the numerical evaluation of irreducible two-loop corrections. 



In Figure ( TTol l we comp are the im p leme ntation of i?had(s) taken from Burkhardt IU32I1 (i?had,i) and our 
parametrization based on Il47llll49ll ll50ll (flh a H tt). As already stated, the deviations are evidently much 
smaller than one might expect and may be considered to be irrelevant here. 

We close this section with a brief discussion of narrow resonances. Narrow resonances are implemented 
replacing the rapidly varying cross section ratio with the parametrization 



9tt 



Rr C s(z) = ^M ros r^ es e 



S(z - M4) 



(E3) 



The integration over z is then carried on analytically leading to the following result for the IR-finite remainder 
(including the irreducible box diagrams) of Eq. ( |87| i: 
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F 2 (Af 2 cs ) + F 3 (Af 2 GS )ln 
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(E4) 



For the numerical evaluation of the contribution due to the narrow resonances, we use the values listed in the 
Burkhardt' s routine 113211 . collected in Table W\ 



APPENDIX F: EVALUATION OF POLYLOGARITHMS 



At several instances, dilogarithms \A<z(z) and trilogarithms ~L\${z) of complex argument are needed. A 
definition of polylogarithms is: 



Li n (z) = S n -i,i{z) — 
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FIG. 16: A comparison of the parametrizations from U3a1 and U4^\] . 



resonance 


Mres [GeV] 


r r e e ; e " [keV] 


oj(782) 


0.7826 


0.66 


<£(1020) 


1.0195 


1.31 


JMiS) 


3.0969 


4.7 


V(2S) 


3.6860 


2.1 


V>(3770) 


3.7699 


0.26 


V>(4040) 


4.0300 


0.75 


^(4160) 


4.1590 


0.77 


^(44 15) 


4.4150 


0.47 


T(1S) 


9.4600 


1.22 


T(2S) 


10.0234 


0.54 


T(3S) 


10.3555 


0.40 


T(4S) 


10.577 


0.24 


T( 10860) 


10.865 


0.31 


T( 11020) 


11.019 


0.13 



TABLE V: Numerical values for the treatment of narrow resonances, taken directly from 113 



They have the special values Li n (0) = and Li n (l) = ((n), where £(s) is the Riemann ^-function, £(2) = 
7r 2 /6,C(3) = 1.2020569031595942854. . . An efficient evaluation transforms the arguments to the region 
where modulus and real part are bound: \z\ < 1 and die(z) < i, using: 



(M-iln 2 (-*)-C(2), 



Li 2 (z) 



-Li 2 (l -z) + C(2) - \n(z) ln(l - z), 



(F2) 
(F3) 
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and 



U 3 (z) = Li 3 Qj-iln 3 (-z)-C(2)ln(-z), 

Lia(z) = -Li 3 (l " ;) - Li 3(l - *) + C(3) + i ln 3 (z) + C(2) ln(*) - i ln 2 (z) ln(l - z) 



(F4) 
(F5) 



Th en, se ries expansions with Bernoulli numbers ensure rapid convergence. For Li2(z) we follow Appendix A 
oflflSltl: 



Li 2 (z) = £ 
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AC; + Di 



■ ln(l - «)] 
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1 o 

ln(l-z)- -In (l-z) + 



4^EC(2J) 



i=i 



(-iy 

2j + l 



ln(l - z) 



1 2j + l 



2?r 



(F6) 



The _Bj are B ernoulli numbers, Bq = 1, etc. Useful series expansions for Li n (z) are given in Eqns. (48) and 
(49) of Il52ll . which we reproduce here for the special case n = 3: 



U 3 (z) = J2 



C 3 (j) 



3=0 

j 



c 3 (j) = £ 



fc=0 



(?+l)l 

f j ) B 3 -kB k 
1 + fc ' 



ln(l - z)] 



j'+i 



(F7) 



(F8) 



with Cs(O) = 1 etc. For Li2(z) and Li 3 (z) we observe typically that n summation terms give an n ± 1 digits 
accuracy. We just mention that we do not allow to evaluate the logarithms and polylogarithms at their cuts 
(negative real axis beginning at z = and positive real axis b eginn ing at z = 1, respectively). For other 
conventions we refer to the c orres ponding remark at p. 19 of Il52ll . Our Fortran code is available as file 
cpolylog . f at the website llll7ll . 

An alternative, efficient algorithm for the evaluation of polylogarithms is described in lll53ll lll57ll . 
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